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IMPINGEMENT OF WATER DROPLETS ON WEDGES AND DOUBLE-WEDGE AIRFOILS 

AT SUPERSONIC SPEEDS ^ 

By John S. SERAriNi 


SUMMARY 

An analytical solution is presented for the equations of motion 
of water droplets impinging on a wedge in a two-dimensional 
supersonic flow field with an attached shock wave. The closed- 
form solution yields analytical expressions for the equation of 
the droplet trajectory, the local rate of impingement and the im- 
pingement velocity at any point on the wedge surface, and the 
total rate of impingement. The analytical expressions are 
utilized in the determination of the impingement of water drop- 
lets on Hue forward surfaces of symmetrical double-wedge aiifoils 
in supersonic flow fields with attached shock waves. 

For a wedge, the results provide information on the efects of 
the droplet size, the free-stream Mach number, the semiapex 
angle, and the pressure altitude. For the douile-wedge airfoil, 
additional calculations provide information on the effect of air- 
foil thickness ratio, chord length, and angle of attack. 

The results for the symmetrical double-wedge airfoil are also 
correlated in terms of the total collection fficiency as a function 
of a relative modified inertia parameter. The results are pre- 
sented for the following ranges of variables: free-stream static 
temperature, 4^0° to 4^0° B; droplet diameter, 2 to 100 microns; 
free-stream Mach number, 1.1 to 2.0; semiapex angle for the 
wedge, 1.14° to 7.97°, and corresponding double-wedge-aitfoil 
tidckness-to-chord ratio, 0.02 to 0 . 14 ; pressure altitude, sea 
level to 30,000 feet; and chord length, 1 to 20 feet. 

INTRODUCTION 

The problem of ice prevehtdoa on aircraft fljniig at sub- 
sonic speeds up to critical Mach numbers has been a subject 
of considerable study and research by the NACA. The recent 
advent of aircraft flying at transonic and supersonic speeds 
has required an extension of these icing studies. That an 
icing problem exists in the transonic and supersonic speed 
range is verifled in reference 1, which shows by an analytical 
investigation with experimental confirmation that diamond 
or symmetrical double-wedge airfoils are subject to possible 
icing at flight Mach numbers as high as 1.4. A similar re- 
sult is expected for other airfoil shapes being considered for 
use at transom’c and supersonic flight speeds. 

In conducting research on the problem of ice prevention 
on aircraft and missiles, regardless of the magnitude of the 
flight speed, it is essential that the impingement of cloud 
droplets on airfoils and other aerodynamic bodies be deter- 
mined either by theoretical calculations or experiment. The 


impingement variables that must be determined are the total 
water catch, the extent of impingement, and the rate of im- 
pingement per unit area of body surface. These variables 
can be determined analytically from calciflations of the cloud- 
droplet trajectories obtained for the various aerodynamic 
bodies. Investigators have reported the results of studies of 
doud-droplet trajectories about right-circiflar cylinders (refs. 
2 to 5) and about airfoils (refs. 6 to 9) immersed in an in- 
compressible fluid. An evaluation of the efliect of compres- 
^sibility on the droplet trajectories about cylinders and airfoils 
up to the critical flight Mach number is presented in refer- 
ence 10. 

At present, little infonnation exists on the impingement 
of droplets on aerodynamic bodies in a supersonic air stream. 
The concentration of past effort on problems of impingement 
on airfoils at subsonic flight speeds and the present lack of 
convenient and rapid means for obtaining the rotational 
flow fields about airfoils at supersonic speeds are possible 
explanations for the scarcity of trajectory calculations for 
the supersonic region. An initial contribution to the solu- 
tion of the over-aU problem of impingement of water droplets 
on aerodynamic bodies at supersonic speeds is given in 
reference 11, which presents an analysis of the water- 
interception characteristics of a wedgein asupersonicflowfield. 

The present report extends the analysis of reference 11 
and fmther presents an extensive study of the impingement 
of water droplets on two-dimensional wedges and double- 
wedge airfoils for supersonic flight speeds that result in 
attached shock waves and constant velocity fields behind the 
shock waves. For the wedge angles and double-wedge- 
airfoil thickness ratios considered herein, the shock-wave- 
attachment Mach number varies from a value slightly 
greater than 1 to about 1.4. The method employed is based 
on ah analytical solution of the equations of motion by means 
of a closed-form integration. The closed-form solution 
yields analytical expressions for the equation of the trajec- 
tories, the local impingement efficiencies, the velocity at any 
pomt on the trajectories, and the total rate of impingement. 
This solution is made possible by the use of an empirical 
relation for the drag coefficients for spheres that gives a good 
approximation of the experimental drag coefficients. 

The results of calculations for the rate, the extent, and the 
distribution of impingement of water droplets on wedges 
and symmetrical double-wedge airfoils are presented herein. 


* Sujwsedes NAOA TN 2971, **Impfngcment of Water Droplets pn Wedges and Diamond AlrfoDs at Supersonic Speeds/^ by John 8. Serafinl, 1953. 
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Tlie ranges of variables included for tbe wedge are: free- 
stream static temperature, 420°, 440°, and 460° R; droplet 
diameter, 2 to 100 microns; free-stream Mack number, 1.1 
to 2.0; tangent of tbe semiapex angle, 0.02 to 0.14; and 
pressure altitude, sea level, 15,000, and 30,000 feet. The 
ranges of variables for the double-wedge airfoil are the same 
as those for the wedge, and the additional variables for the 
symmetrical double-wedge airfoil range from 0.02 to 0.14 for 
the thickness ratio and from 1 to 20 feet for the chord length. 

The work reported herein was performed at the NACA 
Lewis laboratory in the spring and summer of 1952. 

ANALYSIS 

STATEMENT OF PEOBLEM 

The solution of the problem of impingement of water 
droplets on a two-dimensional wedge at supersonic speeds 
with an attached shock wave is not as difficult* as that for 
the impingement on various airfoils at low subsonic speeds. 
For the wedge at supersonic speeds with an attached shock 
wave, the air velocity everywhere ahead of the shock wave 
is constant and equal to the free-stream air velocity Vi 
(fig. 1). The air velocity behind the shock wave Vj is also 
everywhere constant and parallel to the wedge surface. All 
the droplets have the same initial velocity (that of the free-* 
stream air velocity), and their trajectories are exactly co- 
incident with the air streamlines upstream of the shock wave. 
All water droplets of a given size are subjected to identical 


air-velocity fields, which in turn produce identical force 
systems downstream of the shock wave, hrespectivo of the 
point along the shock wave where the droplets cross the 
wave. It follows, therefore, that, for droplets of a given 
size, all the trajectories in a given problem are identical with 
respect to the point where the droplet crosses the shock wave. 

By adopting a frame of reference that moves at the con- 
stant velocity of the air Vj downstream of the shock wave, 
the problem of the droplet motion is reduced to the stUl-air 
problem, defined as the' determination of the motion of a 
droplet that, having an initial velocity, is projected into 
quiescent Hence, relative to the moving frame of refer- 
ence, the initial velocity of the droplet upon crossing the 
shock wave is equal to the vectorial difference of the free- 
stream air velocity Yi and the air velocity Vi doAvnstream of 
the shock wave. Adoption of the frame of reference moving 
with a constant velocity reduces the problem from the solu- 
tion of two simultaneous nonlinear second-order differential 
equations in the fixed coordinate system to the solution of a 
single nonlinear second-order differential equation in the 
moving coordinate system. At any given instant, the drop- 
let displacement relative to the point of intersection wth the 
shock wave in the fixed frame of reference is obtained by 
adding vectoriaUy the droplet displacement ^vlthin the mov- 
ing frame of reference to the displacement of the moving 
frame of reference for the same increment of tune. 


At / = o . ■ At / = / 



Fiqueb 1. — Schematic diagram of water-droplet trajectory for wedge in supersonic flow with attached shock wave. 
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This general method was used in reference 11, where the 
one second -order differential equation representing the drop- 
let motion relative to the air velocity behind the shock wave 
was integrated graphically. However, it is possible to 
obtain a completely analytical solution by means of a closed- 
form integration without resorting to the use of numerical 
integrations or analog computing equipment, if an empirical 
relation is assumed for the drag coefficient as a function of the 
Reynolds number of the droplet relative to the air. It will 
be shown that this closed-form integration of the still-air 
problem when applied to the wedge in supersonic flow with 
attached shock wave yields the equations for the trajectories 
of the water droplets and the droplet velocities at any point 
on the trajectories and makes available relations for the 
rates of total water impingement and the local rates of water 
impingement along the wedge surface. Furthermore, it is 
shown that these equations can also be readily applied to the 
determination of the droplet impingement on a double-wedge 
airfoil in supersonic flow with attached shock wave. 

Three of the usual assumptions made in the previous in- 
vestigations on impingement at subsonic speeds and also 
required for this investigation are (1) the water droplets are 
always spherical and do not change in size, (2) the force of 
gravity on the droplet may be neglected in comparison with 
the drag forces, and (3) the drag of the air on the droplet is 
that of a viscous incompressible fluid. Here it is additionally 
assumed that (4) the two-dimensional supersonic flow field 
about the wedge is frictionless except within the infinitesi- 
mally thin attached shock wave, (5) no condensation shock 
occurs and no change in phase occurs as the water droplets 
traverse the oblique shock wave, and (6) the imbalance of 
the forces on the water droplet from the instant it enters the 
shock wave until it emerges from the shock wave can be neg- 
lected in the calculation of the trajectories. 


EQUATION OF DBOPLET MOTION IN MOYINQ KBFEBENCE FBAME 

The velocity of the droplet in the moving frame of refer- 
ence is 

U=\Va-V,\ (1) 


where Va is tha droplet velocity with respect to the fixed 
frame of reference, and Ta is the air velocity downstream of 
the attached shock wave also with respect to the fixed 
frame of reference (fig. 1). In the frame of reference moving 
with the velocity Vt, the statement of Newton’s law of mo- 
tion for the water droplet becomes 

D=-Oo^ ^ ( 2 ) 


from which 


dU_ _ 3 p2 W p 
dt~~ 8 a 


( 3 ) 


(A complete list of symbols is given in appendix A.) 

Equation (3) is the differential equation of motion of a 
droplet projected with an initial velocity into a region of 
quiescent air (the so-called stfll-air problem). The shock 
wave is considered to be a surface of discontinuity from which 


the droplets emerge with a velocity "Fj. In this case the 
initial velocity of the droplet is 

Ui=\Vi-V,\ (4) 

which is the magnitude of the vector difference of the air- 
velocity vectors upstream and downstream of the attached 
shock wave. As can be shown from consideration of the 
continuity equation and the equation for conservation of 
momentum across the oblique shock wave, the velocity vec- 
tor Ut is normal to the shock wave. At any subsequent 
instant of time, the relative droplet velocity vector U retains 
the same angular orientation to the shock wave and changes 
only in magnitude. 

In reference 11 the solution of equation (3) is obtained by 
numerical integration. The result obtained in this manner 
makes it necessary to use a graphical procedure in determin- 
ing the trajectories and the local rates of impingement. 
However, an analytical solution of equation (2), which elimi- 
nates the graphical procedure, can be obtained if the experi- 
mental values of the drag coefficient Gd are expressed in a 
function involving the Reynolds number Rer- The relation 
is 


where e and m are the empirical constants. This empirical 
relation is a valid approximation in the range of Reynolds 
numbers to which cloud droplets are subjected in trajectory 
calculations. Substitution of the expression for Go (eq. (5)) 
in equation (3) results in the expression - 


dU _£ ^ H 

dt dt^ 2 Pa a* 




( 6 ) 


where the Iqcal relative Reynolds number Rer=2p2UaJpi. 
The displacement of the water droplet in the moving frame 
of reference x is measured from the air streamline that inter- 
sects the shock wave at the point where the water droplet 
enters the air-flow field downstream of the ^ock wave. 
The closed-form integration of the differential equation (6) 
is presented in appendix B. ■ The use of 2/3 for the exponent 
m and 0.158 for the value of the empirical constant e in equa- 
tion (5) yields an empirical curve for the drag coefficient as 
a function of the local Reynolds numbers that approximates 
very well the variation of the experimental values of the 
drag coefficient hi the range of Reynolds numbers from 0.5 
to 500. The value of 2/3 for the exponent m also facihtates 
the closed-form integration of the differential equation of 
motion. In figure 2 a graph of the empirical relation is 
presented, along with the drag-coefficient data of references 
4 and 12. 

The results of the integration are given by the following 
equations: 


6 P2 L 

V(Ee,. l)e-- 1 V(Ee,.r^^£->-t-l)e-^lJ 

( 7 ) 


308505—01 


-7 
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and 

*»=| ^ (Ber, (9) 

where 

<P=t&n-^ (Rer.r^ih-^'^; 

The intermediate steps of integration are _given in appendix B. 

Equations (7) and (8) give, respectively, the displacement 
and the velocity of the droplet at any instant in the moving 
frame of reference. The displacement of the droplet with 
respect to the point where it crosses the shock wave can be 
obtained by a vectorial addition of the displacement x and 
the displacement of the- moving reference frame in the cor- 
responding time interval. The droplet velocity Vi relative 
to the fixed frame of reference must also he obtained by the 
vectorial addition of TJ (eq. (8)) and V^. Equation (9) 
gives the maximum value of x obtained as the time of travel 
in the air-flow field downstream of the shock wave 
approaches infinity. The significance of this quantity wifi, 
bo discussed in subsequent sections. 


EELATIONS EEQDIRED FOE APPUCATION OF CLOSED-FOHM SOUTTION 
TO OBTAIN DEOPLET MOTION AND IMPINGEMENT IN FIXED 
EEFEEENCE FEAME 

Impingement on wedges. — ^For a problem of given aero- 
dynamic conditions, the trajectories of all the water droplets 
of a given size are identical when the points where the droplet 
trajectories intersect the shock wave are superimposed. 
This unique characteristic of the water-droplet trajectories 
about a wedge in supersonic flow with an attached oblique 
shock wave is the result of two constant velocity fields, one 
upstream and one downstream of the shock wave. There- 
fore, only one set of equations for a single trajectory is neces- 
sary to calculate the impingement parameters for a specified 
problem, including a given droplet size. The values of the 
initial relative velocity Z7<, the initial Reynolds number 
Rer,t, and the density ratio pjpj are needed for substitution 
in the closed-form solution of the equations of droplet motion. 
These values can be obtained from information available in 
reference 13 and from the use of simple algebraic and trigo- 
nometric relations for given values of the free-stream static 
temperature U, the droplet diameter d, the free-stream Mach 
number Mi, the angle of surface inclination to the free-stream 
direction <r, and the free-stream static pressm’e pi. These 
relations result in the following expressions for the initial 
relative velocity and the initial relative Reynolds number: 



Fig ORB 2. — Comparison of empirical relation for drag coefiBcient as function of droplet Reynolds number with experimental values. 
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Ut=\Ui\= FiVl+i2'-2Q COS <r= Viw (1 0) 

and 


S=x sec (y-\-a)-{-V 2 




r 


(15) 


A convenient form of the solution for the impingement on 
a wedge or the front half of a double-wedge airfoil is obtained 
if S is defined as the distance to the point of impingement 
measured from the leading edge for a water droplet that 
enters the flow field behind the shock wave at a distance f 
above the leading edge (fig. 1). The irnique relation between 
S and f in a given problem for droplets of the same size is 
quite readily determined by considering the displacement of 
the Avater droplets as the vectorial sum of the displacement 
of the water droplet relative to the moving frame of reference 
and the displacement of the moving reference frame relative 
to a fixed frame of reference (referred to wedge). Since the 
moving reference frame has a velocity equal to the air velocity 
Va, which is constant in magnitude and parallel to the wedge 
surface, only the droplet travel in the moving reference frame 
includes the component of droplet travel representing the 
approach of the water droplet to the wedge surface. For a 
water droplet starting from point A and impinging on the 
wedge surface at point D (fig. 1), the displacement of the 
moving reference frame (=Fji, where i is zero at p oint A) is 
given by the displacement vector AB equal to CD, and the 
droplet motion in the moving frame of reference is given by 
the displacement vector BD equal to AC. Therefore, rela- 
tive to the starting point at A (fig. 1), the displacement of 
the water droplet to the point of impingement at D is 
obviously equal to ^-hBD or to AC-hCD. From figme 1, 
the displacement of the water droplet at the point of impinge- 
ment D, measured from the leading edge at E, is given by 
adding the vector EA to the displacement vector from the 
starting point A; and the displacement of the droplet at D 
referred to the leading edge is 


<S=|ED|=|EA-t-AC-t-CD|=|EC-t-CD|=|EC|-t-|CD|=f-|-{' 

( 12 ) 


where f is the magnitude of the displacement vector EC, 
and i=Vat) is the magnitude of the displacement vector 
CD (the displacement of the moving frame of reference). 

The values of f and f are obtained in terms of x, the distance 
of travel in the moving reference frame, from simple trigo- 
nometric identities involving the various angles shown in 
figure 1: 

i=x sin 6 tan (r+<r) , (13) 

where x is given by equation (7); r by 

!] 

and 

i=x sec (v+(t) (14) 



Substitution of equation (14) into equation (12) for the sur- 
face distance to the point of impingement yields the following 
equation: 


Since a; is a function of r in equations (13) and (15), the 
expressions for f and S, respectively, are functions of r. 
However, since t cannot he eliminated from equations (13) 
and (15), S cannot he obtained explicitly as a function of f. 
Nevertheless, the curves of f, the initial displacement of the 
water-droplet trajectory from the leading edge normal to 
the free-stream direction, against S, the distance to the point 
of impingement of the stated water-droplet trajectory, can 
be obtained by substitution of the same set of values for r 
in the expressions for f and S. 

An analytical expression for the local impingement 
efficiency can be obtained from the preceding expressions 
for f and S. The local impingement efficiency /3 is defined 
by the expression 

dr 

where is the difference in the initial displacements of two 
water droplets having very nearly equal initial displacements, 
and AS is the small increment of wedge surface between the 
points of impingement of the two water droplets. Differen- 
tiating f and S with respect to t and performing the division 
indicated by equation (16) yield the following expression 


for /3: 

• 



Til 7t\ Cb pfg 

(17) 

where 

7ij=sin 6 tan (v+cr) 

(18a) 


7^3=8ec. (r-hcr) 

(18b) 


-rr -i-T t 

7^3 — F 2 -X — Va — 

opa T 

(18c) 


Since P and S both are fimctions of r, the local impingement 
efficiency /3 at any point on the wedge surface is determined 
by using the same value of t in equations (17) and (15). The 
value of that exists as the point of impingement of the 
water droplet on the wedge surface approaches the leading 
edge as a limit (S'— >0) is given by the following: 


/3o=lim (/3)= 

S-iO 


1 


3 fia 1 

Til Til I? Pa Ui 


= Sm ff 


(19) 


The magnitude and direction of the droplet velocity at the 
point of impingement Vd,im (relative to the fixed frame of 
reference) can also he easily obtained at any point on the 
wedge surface as 

F,*. Fi ^ cos (v-hv) (20) 

/Ci„=o-— o-|„=o-— sin“^ sin (r-i-o-)J (21) 

where k is the angle between the free-stream velocity vector 





IMPINQBMENT OP WATER DROPLETS ON WEDGES AND DOUBLE-WEDGE AIRFOILS AT SUPERSONIC SPEEDS 91 


of trajectories of droplets impinging on the surface beyond 
the shoulder or 60 percent of chord of the double-wedge 
airfoil (surfaces BD or CD), where a portion of the trajectories 
is within the expansion zone emanating from the shoulder. 

RESULTS AND DISCUSSION 

From the equations presented in the previous section and 
in appendix B, the impingement of water droplets on a wedge 
in a supersonic flow fleld with an attached shock wave can 
be calculated over a range of free-stream conditions, wedge 
angles, and droplet sizes. As has been indicated previously, 
the impingement characteristies of a double-wedge airfoil 
(at zero angle of attack and also for small values of angle of 
attack) can readily be determined from the impingement on 
wedges for similar conditions. The results for the wedge 
and the double-wedge airfoil are presented and discussed 
separately. A comparison of the total collection efficiency 
and the water impingement rate at zero angle of attack for a 
double-wedge airfoil with those for an NACA 0006-64 
ahfoil is presented in appendix C. 

WEDGE 

local impingement efficiency. — The rate of water impinge- 
ment on a local area of wedge or airfoil surface is proportional 
to a dimensionless term /3, the local impingement efficiency. 
The local rate of water impingement in poxmds per hour 
per square foot is 

TFiS=0.3296FiW/3 

where /3 is deflned in equation (16). The local impingement 
efficiency, when given as a function of the surface distance 
of the wedge, allows the determination of the local rate of 
impingement of water droplets at any point on the surface, 
the total impingement of water droplets on the entire sur- 
face or any given portion of the siuface, and the extent of 


impingement on the surface. The local impingement 
efficiency /3 is related to a point at a given distance S on 
the wedge surface m equations (15), (17), and (19) by the 
dimensionless time variable r, which is common to all three 
expressions. The variation of /3 with S is presented in 
figure 4 for an extensive range of free-stream conditions, 
semiapex angles, and droplet sizes. The value of J3 at the 
leading edge (S-^O) is the sine of the semiapex angle (sin a) ; 
and as S increases, the value of (3 decreases rapidly and ap- 
proaches the value of zero asymptotically as S approaches 
infinity. However, it is to he noted that negligibly small 
values oi p (fi~l percent of po for the wedge) are attained 
at large but finite values of S. 

The curves of jS as a function of S presented in figure 4 
are those of an idealized situation. The assumed two- 
dimensional supersonic flow field about the wedge does not 
account for a stagnation point that must exist at the leading 
edge of the wedge, regai-dless of the sharpness of the leading 
edge. In addition, the leading edges of wedges and double- 
wedge airfoils might be considered to be somewhat roimded 
when subjeeted to considerable magnification. Therefore, 
it is reasonable to assume that very near the leading edge 
(S-^0) the value of p would actually he greater than the 
calculated value of P at the given distance S. However, 
this should have a neghgible effect on the rest of the P curve 
and also on the total impingement on the wedge, since the 
effect of a stagnation point would bo limited to a very 
small region about the leading edge. 

For S approaching very large values, the calculated values 
of p probably differ somewhat from actual values obtained 
in flight, because the analytical solution of the present report 
does not consider the existence of the boundary layer on the 
wedge surface. Since the boundary-layer thickness increases 
with the surface distance along the wedge, droplets that im- 



(a) Efifeot of pressure altitude. Droplet diameter, 20 microns; free-stream Mach number, 1.4; tangent of semiapex angle, 0.02. 
Fiqueb 4. — Local impingement efficiency on wedge as function of distance along surface. Free-stream static temperature, 440° R. 
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piugo at large distances from the leading edge actually would 
have traveled through the boundary layer for some non- 
negUgible interval of time. However, only a very small 
fraction of the total water droplets of a given size impinge 
under this condition, and for large values of S the values of 
P are negligibly small. For example (fig. 4 (a)), such would 
be the case for values of S greater than 8 or 9 feet, where 
/3<0.0002 as compared with /3 q=0.020. 

A preliminary smwey disclosed a neghgible effect of the 
free-stream static temperature on the local impingement 
efficiency as a fimction of the surface distance along the 
wedge (/S against iS). Values of jS and corresponding values 
of S were calculated for free-stream static temperatures of 
420°, 440°, and 460° R, droplet diameter of 20 microns, free- 
stream hlach number of 1.3, tangent of the semiapex angle 
of 0.06, and pressure altitude of 15,000 feet. The values of 
P for the free-stream static temperatures of 420° and 460° R 
are Avithin 1 percent of the values of p at 440° R. Since these 
calculations show that curves of p agaiost S for the three 
values of free-stream static temperature form practically a 
single curve when plotted to the usual scales, no figure is 
presented to illustrate the effect of free-stream static tem- 
perature on impingement. Furthermore, the results included 
herein, which are calculated for a free-stream static temper- 
ature of 440° R, may be used in the range of temperature 
from 420° to 460° R or possibly an even greater range. 

The effect of the free-stream static pressure on the local 
impingement efficiency as a function of the surface distance 
along the wedge is presented in figure 4 (a) for pressme 
altitudes of sea level, 15,000, and 30,000 feet. Increas i ng 
the pressure altitude (decreasing the free-stream static pres- 
sure) increases the values of p at any distance S. For ex- 
ample, at S=1.5 at sea level, P is 0.0052; and at 30,000 feet, 
P is 0.0071. Since Po (the value of p at i5=0) is equal to 
sin <r, whore o- is the semiapex angle, the curves for the three 
pressure altitudes have the same maximum value of p. Also, 
for the various pressure altitudes, the extent of impingement 
along the wedge surface is essentially the same. 

The effect of the semiapex angle v of the wedge on the 
local impingement efficiency as a function of the surface dis- 
tance along the wedge is presented in figure 4 (b) for values 
of tan V from 0.02 to 0.10. Since the values of Po are equal 
to sin <r, increasing the semiapex angle of the wedge results 
in an increase of po. The surface extent of perceptible im- 
pingement (as characterized by 0.01ft)) does not vary as 
the wedge thickness is increased. 

The effect of free-stream Mach number on the p curve is 
presented in figure 4 (c) for free-stream Mach numbers of 
1.2, 1.3, 1.4, 1.5, and 2.0. For the wedge semiapex angle 
presented in the figure (tan o-=0.04), the value of Mi=1.2 is 
close to the shock-wave-attachment Mach number. The 
shock-wave-attachment Mach number is a function of the 
Avodge semiapex angle and is defined as that Mach number 
below which the shock wave is detached from the wedge. 
An increase in the free-stream Mach number Mi results in 
an increased surface extent of perceptible impingement and 
also in an increased value of p at any given distance S (except 



(b) Effect of wedge semiapex angle. Droplet diameter, 20 microns; 
free-stream Mach number, 1.3; pressure altitude, 16,000 feet. 


Figure 4. — Continued. Local impingement efficiency on wedge os 
function of distance along surface. Free-stream static temperature, 
440° R. 

at S=0, where p=Po^8m <r and at where p—>0). 

This increase in the surface extent of perceptible impinge- 
ment is shown in figure 4 (c), in which, for free-stream Mach 
numbers of 1.2, 1.3, 1.4, 1.5, and 2.0, the surface extents 
of perceptible impingement on the wedge (tvhere /3«0.01/3o) 
are 5.35, 6.05, 6.65, 7.20, and 9.4 feet, respectively. 

The effect of the droplet size on /3 is presented in figme 4 (d) 
for droplet diameters of 10, 20, 30, 40, 50, and 100 
microns. The surface extent of impingement and the values 
of p at any given distance S are considerably increased as the 
droplet size is increased. For example, for the semiapex 
angle presented in the figure (tan a-=0.06), at S=3 feet the 
values of p are 0.0000, 0.0033, 0.0106, 0.0182, 0.0244, and 
0.0417 for values of droplet diameter of 10, 20, 30, 40, 50, 
and 100 microns, respectively. The surface extent of per- 
ceptible impingement has values of 1.5, 5.7, 11.2, 18.1, and 
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FiauEE 4. — Concluded. Local impingement efficiency on wedge as function of distance along surface. Free-stream static temperature, 440° R. 
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28.0 feet for droplet diameters of 10, 20, 30, 40, and 60 
microns, respectively. As shown by the preceding dis- 
cussion and also by a comparison of the values of /3 as a 
function of S (fig. 4), varying the droplet diapieter from 20 
to 30 microns or from 30 to 40 microns is of the same order 
of magnitude in its effect on the /3 against iS curve as varying 
the pressure altitude from sea level to 30,000 feet or varying 
the free-stream Mach number from 1.2 to 1.5. 

Total impingement. — ^The effects of the free-stream Mach 
niunber, the semiapex angle of the wedge, the pressm’O alti- 
tude, and the droplet size on the total impingement on a 
wedge siu4ace of infini te extent are given in figmn 5. The 
total impingement is represented by which is the f of 
the droplet having its trajectory tangent to the wedge smface 
(theoretically the tangent trajectory touches the wedge 
surface at a point S-*-<o ). This can be obtained from the 
integration 

or more directly from equation (13) of the analytical solution, 
as 

sin 0 tan (j'-f-ff) 

where x„ is given by equation (9) and is also defined in appen- 
dix A. The value for is obtained from the expression for 
X (eq. (7)) by allowing t-^». Since the droplet that cntem 
the flow field downstream of the shock wave at a distance 
f equal to (fig. 1) theoretically impinges on the wedge 
surface only as the surface distance S approaches infinity, 
only droplets having values of f equal to or less than will 
impinge on the wedge surface of infinite extent. The rate 
of total water catch on one wedge surface in terms of is 
given as 

W„=0.Z2m„ViW 


where W„ is expressed in pounds per hour per foot of span, 
Vi in miles per hour, and w in grams per cubic meter. 
Therefore, is directly proportional to the rate of total 
water catch on the entire wedge surface and is the rate of 
total water catch on one wedge surface per unit of span, 
free-stream velocity, and hquid-water content (in appropri- 
ate imits). 

The variation of with the tangent of the wedge semiapex 
angle tan <r is shown in figme 5 (a) for free-stream Mach 
numbers of 1.2, 1.3, 1.4, and 2.0. As expected from the 
cmwes of /5 as a function of S, the value of increases as 
tan cr increases. However, the rate of increase in with 
respect to tan <r decreases with an increase in tan tr. As con 
be seen from figure 5 (a), increasing the free-stream Mach 
number increases the value of for a constant value of tan <r. 

The variation of with pressure altitude is presented in 
figure 5 (b) for two Mach numbers. In the range of pressm'o 
altitude from sea level to 30,000 feet, the increase of with 


(a) Effect of tangent of semiapex angle. Droplet diameter, 20 microns ; 
pressure altitude, 16,000 feet. 

(b) Effect of pressure altitude. Droplet diameter, 20 microns. 

(c) Effect 6f droplet size. Pressure altitude, 16,000 feet. 

Figubb 6. — Total impingement rate W„ for wedge of infinite extent. 
Free-stream static temperatvue, 440° R; Tr«=0.329Gf» Viw. 
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an increase in pressure altitude is approximately linear. 

The variation of \vith the droplet diameter d in microns is 
show in figure 6 (c). In the range of droplet diameter from 
10 to 100 microns, as a function of d results in a curve 
that is very nearly a straight line when plotted on logarithmic 
paper. This linearitj’^ permits an accurate interpolation of 
fa, when calculations ai-e made for a few droplet diameters 
for a given value of free-stream Mach number, wedge semi- 
apex angle, and pressure altitude. 

Droplet velpoities at impingement. — ^The variation of 
Vd,imiy\ (ratio of droplet impingement velocity to free- 
stream velocity) -with the surface distance along the wedge 
is presented in figure 6 for three cases. These three cases are 
representative of the results when the droplet diameter d is 
20 microns and the pressm-e altitude is 15,000 or 30,000 feet. 

The curves of as a function of S have characteristics 

similar to the curves presenting /3 as a fvmction of S. At 
S=Q, obwously, all the curves have equal to unity; 

and, as S is increased, the value of the velocity ratio rapidly 
decreases and asymptotically approaches VijVi, the ratio of 
the ah velocity downstream of the shock wave to the ah 
velocity upstream of the shock wave. The sohd curve in 
figure 6 illustrates a tjqiical situation for which V 2 JV 1 is very 
nearly unity (tan <r=0.02 and Mi=2.0). The two lower 
curves (tan (r=0.06 and 0.10 at Mi=1.3) are typical for cases 
where a stronger shock wave produced by a larger semiapex 
angle results in decreased values of Vd tm/Vi for large values 
olS. 

SYMMETRICAL DOUBLE-WEDGE AIRFOIL 

As shoiNTi in the ANALYSIS, the local impingement 
efficiency /3 at any point on the forward surfaces of a double- 
wedge ahfoil (surfaces AB and AC in fig. 3) can be obtained 



Scole parameter, 


(a) Efifeot of pressure altitude and chord length. Droplet diameter, 20 microns; free-stream Mach number, 1.4; airfoil thickness ratio, 0.02. 

FiauEB 7. — Total collection efficiency of symmetrical double-wedge airfoils as function of scale parameter. Free-stream static temperature, 

440° R; angle of attack, 0°. 
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Figure 6. — Variation of ratio of impingement velocity of droplet to 
free-stream velocity with distance along surface of wedge. Free- 
stream static temperature, 440° R; droplet diameter, 20 microns. 


dhectly from the results for the /3 against S curves for wedges. 
The local impingement efficiency )3 may be obtained from 
figure 4 or equation (17). 

In general, the results for the impingement on a sym- 
metrical double-wedge airfoil axe presented in this report in 
terms of the total collection efficiency as a function of the 
scale parameter in an attempt to conform with the 
existing hterature on the impingement characteristics of 
airfoils. In the notation of the present report, the total 
collection efficiency as stated in references 7 and 8 
is defined as 

( 22 ) 


3085C5— 50 8 
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where T is the maxiimiin. thickness of the symmetrical 
double-wedge airfoil, andlfajandlfilare the absolute values of 
the initial displacements from the leading edge (in a direction 
normal to free-stream direction) of the droplet trajectories 
chat impinge at the shoulder of the upper and lower sur- 
faces, respectively, of the double-wedge airfoil. For the 
sjunmetrical double-wedge airfoil at zero angle of attack, |f«| 
and |fi| will be equal; at an angle of attack, the tangent of 
which is equal to the thickness ratio, the value of |f„| is equal 
to zero. 

The scale parameter V' is calculated for the double-wedge 
airfoil at supersonic speeds as for other au-foils at subsonic 
speeds. It is defined as 

9c ^ 23 ) 

a Pa 

where c is the chord length of the double-wedge airfoil. The 
results presented for impingement on double-wedge airfoils 
use essentially the same parameters as used for impinge- 
ment on wedges — droplet size, free-stream Mach number, 
double-wedge-airfoil thickness ratio, and pressm-e altitude. 
In addition, angle of attack and chord length are specified. 

Total collection efficiency at zero angle of attack. — ^The 
variation of the total collection efiiciency E„ with the scale 
parameter ^ is presented in figure 7 for a symmetrical 
double-wedge airfoil at free-stream static temperature of 
440° R and zero angle of attack. The effect of pressure 
altitude on the variation of Em ivith V' for a symmetrical 
double-wedge airfoil of 0.02 thickness ratio is given in figure 
7 (a) for droplet diameter of 20 smicrons and free-stream Mach 
number of 1.4. The pressm-e altitudes presented in the 
figiu*e are sea level, 15,000, and 30,000 feet. The lines for 
constant values of chord length from 1 to 20 feet are also 
included in the figme. For the double-wedge airfoil sub- 
jected to a constant-velocity supersonic flow field, the total 
collection efficiency E„ increases slightly as the pressme 
altitude is increased when the chord length and the other 
variables are held constant. Considering E„ as a function 
of scale parameter ^ where the droplet free-stream JReynolds 
number Rei is held constant yields results s imil ar to those 
for airfoils having rounded leading edges at subsonic speeds. 
With Rci held constant, the value of E„ decreases as ^ in- 
creases in the manner indicated in figme 7 (a). As pre- 
viously mentioned, the total collection efficiency and the 
impingement rate at zero angle of attack for a symmetrical 
double-wedge airfoil are compared in appendix C with 
those for an NACA 0006-64 airfoil. 

The effect of airfoil thickness ratio on the variation of E„ 
vdth ^ is presented in figure 7 (b) for thickness ratios from 
0.02 to 0.14, droplet diameter of 20 microns, free-stream 
Mach munber of 1.4, and pressme altitude of 15,000 feet. 
The droplet free-stream Reynolds number is maintained 
at a value of 453 for all the curves. The effect of increasing 
the airfoil thickness ratio is to decrease the total collection 
efficiency. The rate of decrease in E„ with an increase in 
airfoil thickness ratio becomes somewhat smaller as the 
thickness ratio increases. ' 

The effect of free-stream Mach munber on the variation 
of total collection efficiencgr E„ with scale parameter ^ is 


shown in figme 7 (c) for free-stream Mach munbers of 1.1, 

1.2, 1.3, 1.4, 1.5, and 2.0, droplet diameter of 20 microns, 
airfoil thickness ratio of 0.02, and pressme altitude of 16,000 
feet. The results show that the total collection efficiency 
increases as the free-stream Mach number increases. How- 
ever, the increase in the total collection efficiency from the 
free-stream Mach munber of 1.1 to 1.2 is considerably 
greater than the increase in efficiency from a Mach munber of 
1.2 to 1.3 and from 1.3 to 1.4, and so forth. Calculations 
for thickness ratios of 0.04 and 0.06 indicate the same trend. 
For the 0.02 thickness ratio, the lowest Mach munber pre- 
sented in figme 7 (c) (M=l.l) is quite close to the limiting 
Mach number for shock-wave attachment. Therefore, the 
rate of decrease of total collection efficiency Avith decreases 
in free-stream Mach munber increases as the Mach number 
approaches the shock-wave-attachment Mach munber as a 
limi t. 

The effects of droplet size and chord length on the varia- 
tion of Em with \(f are shown m figmes 7 (d) and (e) . Figure 
7 (d) presents curves of E„ against ^ for a symmetrical double- 
wedge airfoil of 0.02 thickness ratio at a free-stream Mach 
number of 1.4 and a pressme altitude of 16,000 feet. The 
cmwes axe for constant values of droplet diameter (d= 
10, 20, 30, 40, and 60 microns) as well as for constant values 
of chord length (c=l, 2, 4, 8, and 20 ft). Increasing the 
droplet size greatly increases the total collection efficiency. 
For example, at c=8 feet, the values of Em are 0.096, 0.310, 
0.496, 0.626, and 0.711 at droplet diameters of 10, 20, 30, 
40, and 60 microns, respectivdy. The rate of increase in 
the total collection efficiency as the droplet diameter in- 
creases is 1^ for the larger droplet sizes. This effect can 
also be observed in figme 7 (e), in which cmves of E„ as a 
function of ^ are presented for a symmetrical double-wedge 
airfoil of 0.06 thickness ratio for di-oplet diameters of 2, 10, 
20, 30, 40, 60, and 100 microns, free-stream Mach number of 

1.3, and pressme altitude of 15,000 feet. 

A comparison of figmes 7 (b) and (c) (same droplet size 
and pressme altitude) shows that the effect on the total 
collection efficiency E„ of an increase in the free-stream 
Mach munber from 1.1 to 2.0 is, in general, of the same order 
of magnitude as a decrease from 0.14 to 0.02 thickness ratio. 
For example, for ^ of 1790 (c=8 ft) in figme 7 (b), the values * 
of Em decreased from 0.310 to 0.165 for an increase from a 
0.02 to a 0.14 thickness ratio, respectively, a decrease of 
0.145 in the value of E„. For ^|/ of 1790 (c==8 ft) in figure 
7 (c), the values of Em increased from 0.230 to 0.400 for an 
increase in the free-stream Mach nrunber from 1.1 to 2.0, 
respectively, an increase of 0.170 in the value of E„. 

Another comparison of figmes 7 (c) and (d) (same airfoil 
thickness ratio and pressme altitude) shows that the effect 
on the total collection efficiency E„ of an increase in the 
droplet diameter from 10 to 50 microns is much greater than 
an increase in the free-stream Mach number from 1.1 to 2.0. 
For example, for c=8 feet in figme 7 (d), the value of E„ 
increased from 0.095 to 0.710 for an increase in the droplet 
diameter from 10 to 50 microns, an increase of 0.616 in the 
value of Em- As stated previously, for c=8 feet (^=1790) 
in figme 7 (c), the increase in the value of E„ is 0.170 for a 
corresponding increase in the value of the free-stream 
Mach number from 1.1 to 2.0. For a constant value of 
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(b) Eflfecfc of thickness ratio. Droplet diameter, 20 microns; free-etream Mach number, 1.4; pressure altitude, 15,000 feet; free-stream droplet 

Reynolds number, 463. 

(c) Effect of free-stream Mach number. Droplet diameter, 20 microns; airfoil thickness ratio, 0.02; pressure altitude, 15,000 feet. 

Fiodrb 7. — Continued. Total collection eflSciency of symmetrical double-wedge airfoils as function of scale parameter. Free-stream static 

temperature, 440° R; angle of attack, 0°. 
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(d) Efifeot of droplet size and chord length. Free-stream Mach number, 1.4; airfoil thickness ratio, 0.02; pressure altitude, 16,000 feet. 

(e) Effect of droplet size. Free-stream Mach number, 1.3; airfoil thickness ratio, 0.06; pressure altitude, 16,000 feet. 

Figubb 7. — Concluded. Total collection efficiency of symmetrical double-wedge airfoils as function of scale psirameter. Free-stream static 

temperature, 440° R; angle of attack, 0°. 
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chord length, varying the pressure altitude has a relatively 
small effect on the total collection efficiency (fig. 7 (a)). 

Total collection efficiency as function of angle of attack. — 
The previous discussion is concerned with the total col- 
lection efficiencies for the double-wedge airfoil at zero 
angle of attack only. The effect of angle of attack on the 
total collection efficiency is presented in figure 8 for airfoil 
thickness ratios from 0.01 to 0.08, droplet diameter of 20 
microns, free-stream Mach number of 1.4, pressure altitude 
of 15,000 feet, and chord lengths of 1, 2, 4, and 12 feet. 
The range of angle of attack is from zero to tan“’(T/c). 

The total collection efficiency decreases shghtly as the 
angle of attack increases. The slope of the curve of E„ 
against tan a is zero at a=0, by virtue of the symmetry of 
the double-wedge airfoil. As the angle of attack increases, 
the slope of the E„ against tan a curve becomes negative, 
the rate of change becoming more pronoimced as the angle 
of attack increases. The effect of increasing the chord length 
of the double-wedge airfoil is to decrease shghtly the rate of 
decrease of En with respect to tan a. A significant point 
that is well-illustrated in figure 8 is that, for a constant value 
of chord length, there apparently exists an envelope of the 
family of E^ against tan a curves that have the airfoil thick- 
ness ratio Tjc as the parameter. This envelope curve 
presents, in terms of the angle of attack, the maximum total 
collection efficiency that can be obtained for a symmetrical 
double-wedge airfoil of any thickness ratio, where the 
droplet size, free-stream Mach number, pressure altitude, 
and chord length are considered to bo constant. 

For the double-wedge airfoil at a free-stream Mach 
.number greater than the shock-wave-attachment Mach 
number, the decrease in total collection efficiency with an 
increase in angle of attack is opposed to the trend experienced 
by rounded-leading-edgo airfoils at subsonic Mach numbers 
(irre.spective of the symmetry of the airfoil). For the latter 
typo of airfoil at subsonic speeds, the increase in total 
collection efficiency (as defined herein) with increasing angle 
of attack is accounted for by the greatly increased impinge- 
ment that occurs on the lower surfaces of these airfoils at 
angle of attack. The reduction of impingement occurring 
on the upper surfaces of these airfofis at angle of attack is 
more than balanced by the increased impingement occurring 
on the lower surfaces. On the other hand, for a given double- 
wedge airfoil at supersonic speeds with an attached shock 
wave, the rate of increase \vith angle of attack of the impinge- 
ment on the lower surface is less in magnitude than the rate 
of decrease \vith angle of attack of the impingement on the 
upper surface. This general trend for the double-wedge 
airfoil can be explained with the aid of figme 5 (a), in which 
the effect of tan a (tangent of the semiapex angle) on the 
total impingement rate is presented for a wedge of infinite 
extent. 

The forward upper surface of a symmetrical double-wedge 
airfoil effectively acts as the finite portion of one surface of an 
, i nfin ite wedge that is decreasing its semiapex angle (effect- 
ively the thickness of wedge) as the angle of attack of the 
double-wedge airfoil increases. The forward lower surface 
of this type airfoil effectively acts as the finite portion of one 
surface of an infinite wedge that is increasing its semiapex 



FiauRB 8. — Variation of total collection efficiency of symmetrical 
double-wedge airfoils with angle of attack. Free-stream static 
temperature, 440® R; droplet diameter, 20 microns; free-stream ^lach 
number, 1.4; pressure altitude, 16,000 feet. 

angle as the angle of attack increases. In figure 5 (a) the 
increase in fa, which is exactly proportional to the total 
impingement rate, becomes smaller with an increase in tan <r 
(i.e., (^^^/^(tan tr)*<0) for all possible values of the semiapex 
angle a. For example, in figure 5 (a) for the curve of Mi= 1.4, 
the impingement on a symmetrical double-wedge airfoil of 
0.08 tbickness ratio and large chord extent at zero angle of 
attack would be very nearly proportional to twice the value 
of at tan <r of 0.08 (2f„=0.146 ft). The impingement for 
the same airfoil under the same conditions and at an angle of 
attack of 2.291° (tan a=0.04) would be very nearly pro- 
portional to the sum of at tan v of 0.04 and 0.12 
G-«.» = 0.046 -1- 0.0906 = 0. 1365 ft) . 

Total collection efficiency as function of relative modified 
inertia parameter. — ^For air-flow fields that contain a shock 
wave, such as those considered herein, the water droplets 
upon crossing the shock wave suddenly have a velocity 
relative to the air. For the double-wedge airfoil with an 
attached shock wave, this initial relative velocity is the same 
for all droplets entering the air-flow field downstream of the 
' shock wave. As discussed in previous sections, this common 
initial velocity Ui may be considered the initial velocity of 
projection of a droplet in a reference frame having no air 
motion in it and moving relative to the fixed frame of refer- 
ence at the constant velocity of the air downstream of the 
shock wave. Therefore, it is of interest to define a set of 
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inertia parameters, F and F^, based on the motion of tbe 
droplet in this moving reference frame. The relative 
inertia parameter based on the relative velocity Ut is defined 
as 

p U i Pto(F TJ i $ 

^ yi~9 /xic ~ c 

where K is the inertia parameter in the fixed frame of refer- 
ence, defined as 

K=- (24) 

9 piC 

and is that value of the maximum distance of travel 
obtained when Stokes’ law is assumed for the drag force on 
the droplet. 

The relative modified inertia parameter, also based on tbe 
relative velocity Ut, is defined as 



where x„ is defined in equation (9), and experimental values 
are used for the drag force on the droplet. 

For the problem of water-droplet impingement on airfoils 
having rounded leading edges at subsonic speeds, it is im- 
possible to define a similar relative modified inertia param- 
eter. Also, for any airfoil at supersonic speeds, it is impos- 
sible to define a relative modified inertia parameter, unless 
the shock wave from the leading edge is of constant strength 


in the vicinity of the airfoil, as in the case of the double- 
wedge airfoils with attached shock waves. If the sti'ongth 
of the shock wave is not constant, as for the case of a bow 
shock wave preceding a blunt obstacle, then the value of Fo 
is not constant for all the trajectories of a given droplet size. 
Thus, the problem of impingement on wedges and double- 
wedge airfoils at supersonic speeds with attached shock 
waves is unique in that a relative modified inertia parameter 
can be defined and used in correlating the total collection 
efficiencies. 

The correlation of the total collection efficiency with the 
relative modified inertia parameter Fo is shoivn in figure 9 (a) 
for droplet diameters varying from 2 to 100 microns and 
for pressure altitudes of sea level, 16,000, and 30,000 feet. 
For the symmetrical double-wedge airfoil of 0.06 thickness 
ratio at Mi=1.3 and zero angle of attack, the droplet free- 
stream Keynolds number Rei varies from 42.1 to 2104. The 
values of for a given thickness ratio, Mach number, and 
angle of attack generally form the basis for a single curve 
wiSfi a small amount of scatter existing in the higher range 
of value of Fo. In the lower range of Fo, for all values of Ret 
the plotted points have neghgible scatter. The small scatter 
observed is possibly due to the existence of a very slight 
trend of the curves of F« against Fo with Rei for the higher 
values of Fo. 

The effect of increasing the thickness ratio is presented in 
figure 9 (b) for symmetrical double-wedge airfoils of 0.02, 
0.06, and 0.12 thickness ratio. Increasing the thickness 
ratio displaces the E„ against Fo curve toward larger values 



(a) Correlation of pressure altitude and droplet size. Free-stream Mach number, 1.3; airfoil thickness ratio, 0.06. 

Figubb 9. — Total collection efficiency of symmetrical double-wedge airfoils as function of relative modified inertia parameter. Free-stream static 

temperature, 440° R; angle of attack, 0°. 
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(c) EflFect of free-etream Mach number. Droplet diameter, 20 microns; airfoil thickness ratio, 0.02; pressure altitude, 15,000 feet. 

Figure 9. — Concluded. Total collection eflaoiency of symmetrical double-wedge airfoils as function of relative modified inertia parameter. Free- 

stream static temperature, 440“ R; angle of attack, 0°. 
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of jPo- Changing the thickness ratio of the double-wedge 
airfoil does not alter the shape of the curve itseK. 

The effect of the free-stream Mach number on the variation 
of E„ with Fo is presented in figure 9 (c) for the symmetrical 
double-wedge airfoil of 0.02 thickness ratio for droplet diam- 
eter of 20 microns, free-stream Mach numbers of 1.1, 1.2, 
1.5, and 2.0, and pressure altitude of 15,000 feet. Increasing 
the Mach number displaces the entire curve of E„ against 
Fq toward smaller values of Fq- As the value of Mi increases, 
the rate of displacement of the curve with increasing Mi 
becomes smaller. The displacement of the cmwe obtained 
by increasing the Mach number from 1.1 to 1.2 is more than 
that obtained by increasing it from 1.2 to 1.5 and from 1.5 
to 2.0. 

The relative inertia parameter F and the relative modified 
inertia parameter Fo of the moving frame of reference cor- 
respond, respectively, to the inertia parameter K and a 
modified inertia parameter Kq of the fixed frame of reference. 
The modified inertia parameter is defined in reference 14 as 

.s-«=^a: (25) 

where X is the maximum distance of travel of a droplet 
projected into still air with the free-stream velocity Vi. 
The term X, is the value of the maximum distance of travel X 
when Stokes’ law is assumed for the drag force on the droplet. 
The total collection efficiency of a symmetrical double-wedge 
airfoil is presented in appendix D as a function of the modified 
inertia parameter E^. 

SUMMARY OF RESULTS 

This report presents an analysis of the problem of im- 
pingement of water droplets on a wedge and a double-wedge 
airfoil at supersonic speeds with attached shock waves. 
When a suitable empirical relation is used for the drag co- 
efficient of a sphere, the analysis allows a closed-form inte- 
gration of the equations of motion for the water droplets. 
The integration results in analytical expressions for the 
equation of the trajectories, the droplet velocity at any point 
on the trajectories, the local impingement efficiencies, and 
the total rate of impingement. The results of the calcula- 
tions of rate, extent, and distribution of impingement of 
water droplets on wedges and symmetrical double-wedge 
airfoils are summarized briefly as follows (the Mach number 
referred to is the free-stream Mach number, which is greater 
than the attachment Macti number for the wedge or the 
double-wedge airfoil) : 

1. At a given value of Mach number, droplet size, and 
pressure altitude, the local impingement efficiency as a 
function of the dimensional surface distance is the same for 
both the wedge and the symmetrical double-wedge airfoil 
at zero angle of attack, provided the tangent of the semiapex 
angle of the wedge is equal to the thickness ratio of the sym- 
metrical double-wedge airfoil. 


2. For any Mach munber, pressure altitudei, and droplet 
diameter, the value of j3o is equal to the sine of the semiapox 
angle for wedge or symmetrical double-wedge auloil. [flo 
is the value (maximum) of local impmgement efficiency 
as distance from leading edge to point of impingement 
approaches zero.] 

3. The- effect of the free-stream static temperature on the 
local impmgement effieiency and total collection efficiency 
is negligible for temperatures from 420° to 460° R. 

4. At constant values of Mach number, droplet size, and 
semiapex angle of the wedge or corresponding thickness 
ratio of the symmetrical double-wedge airfoil, an increase in 
the pressure altitude increases slightly the local unpingement 
efficiencies and total collection rates on wedges and S 3 an- 
metrical double-wedge airfoils, but has a negligible effect 
on the surface extent of perceptible impingement. 

5. At constant values of Mach number, droplet size, and 
pressure altitude, increasing the thickness ratio of the sym- 
metrical double-wedge airfoil or corresponding semiapex 
angle of the wedge increases the local unpingement efficiency, 
has a negligible effect on the surface extent of perceptible 
impingement, and decreases the total collection efficiency 
of the symmetrical double-wedge airfoil. 

6. At constant values of droplet size, pressure altitude, 
and semiapex angle of the wedge or thickness mtio of the 
symmetrical double-wedge airfoil, an increase in Mach 
number increases both the surface extent of impingement 
and the value of the local impingement efficiency. 

7. At constant values of pressure altitude, semiapox 
angle of the wedge or thickness ratio of the symmetrical 
double-wedge airfoil, and Mach munber, an increase in 
the droplet size increases considerably the surface extent of 
perceptible impingement, the local unpingement efficiency, 
and the total impingement rate. 

8. The variation of total collection efficiency of the sym- 
metrical double-wedge airfoil at zero angle of attack as a 
function of the scale parameter for constant values of the 
droplet free-stream Reynolds number is similar in form to 
that for subsonic airfoils. 

9. The total collection efficiency of the symmetrical 
double-wedge airfoil decreases slightly as the angle of attack 
increases. 

10. For a symmetrical double-wedge airfoil of a given 
thickness ratio and Mach number, the values of total col- 
lection efficiency for a wide range of values of droplet free- 
stream Reynolds number comprise a single cmwe when 
plotted against the relative modifled inertia parameter. The 
effect of increasing the thickness ratio or decreasing the 
Mach munber is to displace the entire curve in the direction 
of larger values of relative modified inertia parameter. 


Lewis Flight Propulsion Laboratory 
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Cleveland, Ohio, ApT^ 21, 195S 



IMPINGEMENT OF "WATEE DBOPLETS ON WEDGES AND DOUBLE-WEDGE AIKFOILS AT SUPEBSONIC SPEEDS 103 


APPENDIX A 


The 


a 

Co 

c 

D 

d 


F 

Fo 

g 

K 

Ko 

M 

m 

til, tii, 

V 

n 

Re, 

Rci 

S 


following symbols are used, in this report: 
droplet radius, ft 
drag coefficient, dimensionless 
airfoil chord length, ft 
drag force on spherical water droplet, lb 
droplet diameter, microns 
total collection efficiency (defined by eq. 
dimensionless 

relative inertia parameter, ^ 

9 MiC 


SYMBOLS 

7 
€ 


t 

u 

u 


F 

IF 

w 


)3 


(22)), 


dimen- 


sionless 

relative modified inertia parameter. 


^m,t 

,2 


c 


tl3 


acceleration due to gravity, 32.2 ft/sec‘ 

X, 


inertia parameter, 

’ 9 MiC 


dimensionless 


modified inertia parameter, ^K=—i dimensionless 

Aj C 

Mach number 

empirical constant (used in eq. (5)), 2/3 
constants of flow field (defined by eqs. (18a), (18b), 
and (18c), respectively) 
static pressure, Ib/sq ft 
gas constant, 53.3 ft-lb/(lb) (°F) 
droplet Reynolds number relative to air behind 
shock wave, 2o.pj(7/^a, dimensionless 
free-stream droplet Reynolds munber, 2apiFi//i,, 
dimensionless 

distance to point of impingement measiured from 
leading edge for water droplet that enters flow 
field behind shock wave at distance above 
leading edge (eq. (15) and fig. 1), ft 
maximum thickness of airfoil, ft 
time, sec 

free-stream static temperature, °R 
magnitude of droplet velocity relative to air 
velocity downstream of shock wave, \ Vi— Fj|, ]^s 
velocity, fps or mph 
rate of water catch, lb/ (hr) (ft span) 
local rate of water catch, lb/(hx)(sq ft sm^face) 
hquid-water content of air, g/cu m 

J* U dt, displacement of droplet in moving frame 

of reference (relative to air stream); where t=0 
the instant water droplet intercepts shock wave 
maximum value of displacement x as t-> ~ as limit 
(maximum distance of travel of droplet projected 
into still air with relative velocity f7<), ft 
value of maximum distance of travel x„ obtained 
by assuming Stokes’ law for drag force on 
droplet, ft 

angle of attack of airfoil, deg 
local impingement efficiency, dtIdS, dimensionless 
maxi m u m value of as S-^; /3o=sin cr, dimensions 
less 


ir«i, ir.l 


X, 


ratio of specific heats (1.4 for air) 
empirical constant (used in eq. (5)), 0.158, dimen- 
sionless 

initial displacement of droplet trajectories from 
leading edge in direction normal to free-stream 
direction (eq. (13)) 

maximum value of initial droplet displacement f 
obtained when trajectory is tangent to wedge 
surface (theoretically as S—>co) 
absolute values of initial displacement from leading 
edge (in direction normal to free-stream direc- 
tion) of droplet trajectories that impinge at 
shoulder of upper and lower surfaces, respec- 
tively, of double-wedge airfoil 
distance along shock wave measured from wedge 
apex to point where droplet trajectory intercepts 
shock wave 
shock-wave angle 

angle between free-stream velocity vector Vi and 
droplet velocity vector 

maximum distance of travel of droplet projected 
into stOl air with free-stream velocity Fj, ft 
value of distance of travel X obtained by assuming 
Stokes' law for drag force on water droplet, ft 
dynamic viscosity based on static temperature, 
(lb)(sec)/sq ft 

angle between free-stream velocity vector and 
initial relative velocity vector Vt (defined by eq. 
(13a)) 

components of droplet displacement referred to 
wedge surface (defined by eq. (12)), ft 
mass density of air, slugs/cu ft 
water-droplet mass density, 1.9398 slugs/cu ft 
semiapex angle 

angle between droplet velocity vector Vt and 
velocity vector Fj 

dimensionless time variable, {Ziiih>aP?)t 

phase angle, tan“HRer.r'^®e~'^^), (eq. (9)) 

scale parameter, (9c/a) (pjpa) 
ratio of air velocity downstream of shock wave to 
free-stream velocity VijVi 

o) ratio of initial droplet relative velocity to free- 

stream velocity UJVi 
Subscripts: • 

d 
i 

im 

1 

m 
u 
1 

2 


p 

Pw 

(X 


<P 

Q 


droplet 

initial (at shock wave) 
impingement 
lower 
maximum 
upper 
free stream 

downstream of shock wave (in fixed frame of 
reference) 

Barred symbols denote vectorial quantities. 
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APPENDIX B 


CLOSED-FORM INTEGRATION OF EQUATION OF MOTION (RELATIVE TO AIR FLOW DOWNSTREAM OF SHOCK WAVE) 
The various steps necessary to the closed-form integration 


of equation (6), which is the differential equation of motion 
of the water droplet in the moving reference frame, are given 
herein. The differential equation of motion (6) can be re- 
written as 

dt A 
where, for convenience. 


* (i+6Arc7"o 

(Bl) 

■1 a 

Rer=AiUa 

(B2) 


(B3) 


•^2 Q 

0 Ptt 

(B4) 


Algebraic simplification and rewriting result in the form 

As dU 

— m 

where 


'U{l+A,U^) 


and 


A =:4i 
Ai 24 

A^=eA7a” 


(B5) 

(B6) 

(B7) 


Equation (B5) is not readily int^ 3 :able in its present form. 
By letting 

p=J7” 

and 

U” 

"UT 

equation (B5) becomes 


Pi 


(B8) 

(B9) 


dt= 


_ As dp' Aj A^p^ dp' 
m p' ' m l->r A^ip' 


(BIO) 


This form of the equation is readily int^rated, and upon 
resubstituting the relations (B8) and (B9), there is obtained 
the following expression for the velocity of the droplet as a 
function of the time t: 

U= Z7,|exp [((-B.) (Bll) 

where Bi is the constant of integration. 
d'X’ * 

Since equation (Bll) can be rewritten in integral 

form as 

J dx=j Ut |exp [((-B.) ^-Z7TA,p'”d(-l-B, (Bl2) 

where Bj is the second integration constant. Equation 
(Bl2) requires simplification before a closed-form integration 
can be performed. Consider the following substitutions; 


2/=^( (B13) 

As=^t7,(Z7rA,)->/"=^(A0-‘'’" (B14) 

A,=(P7Ad-^ exp [-Bim/A,] (B16) 

By using the substitutions given by (B13), (B14), and (B16), 
it is possible to write equation (B12) as 


J* dx=J Ai(Aie''—l)-^'"dy+Bi 


(B16) 


Before further steps can be taken in the closed-form integra- 
tion of the equation of motion of the droplet, the value of 
the empirical constant m must be determined. It is noted 
that in approximating the curve of the drag coefficient as a 
function of the local Reynolds munber by a relation of tlie 
form given by 

<?i>=^(H-el20 (6) 

it is possible to consider that the value of m is 2/3 and the 
value of € is 0.158. That the approximation of the experi- 
mental cmwe by the empirical relation is very good in the 
range of Reynolds numbers from about 0.6 to 600 can be 
seen by referring to figure 2, which presents a graph of the 
empirical relation along with the drag-coefficient data of 
references 4 and 12. 

The use of ?n=2/3 in equation (Bl6) along -with the sub- 
stitution 

e«'=g (B17) 

and the use of formulas of integration given on pages 16 and 
17 of reference 16 allow equation (B16) to be integrated as 
shown in the following steps: 

Ajdq 


jdz—j 




Ajdq 


q{Atq- 
x—Bi 


.l)i/s 

1 


Aidq 

.l)3/s- 


-B, 


2Ac -^Agq-A 


r Agd^ 

“J (Aeff-l 

-tan“‘ VAeg— 1 (B18) 


If the original time or independent variable t is reintroduced 
into equation (Bl8) and the substitutions are made for Aj 
and A*, the equation has the form 




(A 4 


3A, 


(B19) 
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Substituting for A 3 and (except in the exponent of e) in 
equation (Bl9) results in 


X/ n 3/2 Pl\_ 


a 

x=— 

a V pj 


3 


(B20) 

Equation (B20) is the integrated equation with undeter- 
mined integration constants for the motion of the water 
droplets relatwe to the air 'velocity behind the shock wave. 
The integration constants are determined from the boimdary 
conditions, which are 


x —0 ) 


The substitution of the boundary conditions, and thus the 
determination of the integration constants Bi and B 3 , 
results in the final form of the integrated equation of motion 
for the water droplets relative to the air velocity down- 


stream of the shock wave as follows (in the frame of refer- 
ence moving with the constant velocity Va) : 






( 7 ) 


The final form of the corresponding equation for the relative 
velocity of the water droplet as a function of the dimension- 
less time variable is obtained from equation (Bll) as 


Ui 




( 8 ) 


It can be seen from equations (7) and (8) that as t approaches 
infini ty the value of U approaches zero and that a limit 
exists for the value of a: as i approaches infini ty This 
limiting value of x is 


where 






( 9 ) 


APPENDIX C 

COMPARISON OF COLLECTION EFFICIENCY AND IMPINGEMENT RATE FOE SYMMETRICAL DOUBLE.-WEDGE AIRFOIL 

WITH THOSE FOR NACA 0006-64 AIRFOIL 


A comparison of the total collection efficiency as a func- 
tion of a modified inertia parameter IQ) for a symmetrical 
double-wedge airfoil at supersonic speeds (attached shock 
wave) with that for the NACA 0006-64 airfoil (ref. 9) at 
free-stream !Mach numbers less than critical is presented in 
figure 10. Both airfoils have 0.06 thickness ratios and are 
at zero angle of attack. It is important to keep in mind 
that figure 10 does not allo'w a comparison of the two air- 
foils at the same Mach munber, and such a comparison can- 



Fiqure 10. — CompariBon of total collection efficiency as function of 
modified inertia parameter for symmetrical double-wedge airfoil at 
supersonic speeds with attached shock wave with that for NACA 
0006-04 airfoil at free-stream Mach number less than critical. Air- 
foil thickness ratio, 0.00; angle of attack, 0°. 


not be made, because the analysis for the NACA 0006-64 
airfoil is not valid above the critical Mach number and that 
for the symmetrical double-wedge airfoil is not valid below 
the shock-wave-attachment Mach number. It is, however, 
of value to consider a comparison of the two airfoils -with 
each operating tvithin its appropriate speed range. In figm-e 
10 the rate of increase in E„ with respect to .Ei is greater 
for the symmetrical double-wedge airfoil than for the NACA 
0006—64. For a constant droplet size, pressure altitude and 
temperature, and free-stream Mach number, resulting in a 
constant value of X, varies inversely as the chord length 
iKo=\lc). In general, therefore, the rate of decrease in E„ 
■with increasing chord length is greater in magnitude for the 
symmetrical double-wedge airfoil than for the NACA 0006-64 
airfoil. 

A comparison of the rate of total water catch per uni t 
span for the symmetrical double-wedge airfoil 'with that for 
the NACA 0006—64 airfoil can be obtained if values (neces- 
sarily different) for the free-stream Mach numbers for the 
two airfoils are assigned. Assume that both airfoils are of 
the same thickness ratio and chord length and that they 
encounter identical icing -conditions of droplet diameter (20 
microns), pressure altitude and temperature (15,000 ft and 
440° R), and hquid-water content (0.5 g/cu m), but have 
free-stream Mach numbers of 1.5 and 0.75, respectively. 
For the static temperature assumed, the Mach numbers of 
1.5 and 0.75 correspond to speeds of 1051 and 626 miles 
per hour, respectively. Therefore, the magnitudes of the 
inertia parameter K and the free-stream droplet Reynolds 
number Rei for the double-wedge airfoil are twice as great 
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(a) (Correlation of pressure altitude and droplet size. Free-stream Mach number, 1.3; airfoil thickness ratio, 0.06. 

(b) Effect of thickness ratio. ' Droplet diameter, 20 microns; free-stream Mach number, 1.4; pressure altitude, 16,000 feet. 

(c) Effect of free-stream Mach number. Droplet diameter, 20 microns; airfoil thickness ratio, 0.02; pressure altitude, 16,000 feet. 

Figubb 11. — Total collection efficiency of symmetrical double-wedge airfoils as function of modified inertia parameter. Free-stream static 

temperature, 440° R; angle of attack, 0°. 
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as for the NACA 0006-64 airfoil. Varying the chord length 
from 1 to 20 feet produces a change in the value of from 
0.386 to 0.0193 for the double-wedge airfoil and from 0.266 
to 0.0130 for the NACA 0006-64 airfoil. This variation in 
Kf, for both airfoils results in values of E„ of the same order 
of magnitude. For the given icing conditions, the following 
table Usts for chord lengths of 4 and 20 feet the various per- 
tinent parameters and variables, including the rate of total 
water catch on the airfoil per unit span for both airfoils: 


Symmetrical doable-wedge airfon at Mi—lJi (Vi=10il mph) 

Chord 
length, 
c, ft 

Free-stream 

droplet 

Eoynolds 

number, 

Iffi 

(a) 

Inertia 

parameter, 

K 

Distance 

ratio. 

Modified 

inertia 

parameter, 

J?b 

(b) 

Total col- 
lection 
effldenoy, 

JEm 

Rate of 
water 
catch, Tl^, 
lb/(hr) 
(nnJtspan) 

4 

485 

0.&42 

0.178 

0.0965 

0.432 

17.9 

20 

485 

.108 

.178 

.0193 

.lU 

23.0 

. 

NACA 0006-64 alrioa at A/i-0.7« (V,= 

626 mpb) 


4 

243 

0.271 

0.246 

0.0666 

a330 

6.85 

20 

243 

.0M2 

.246 

.0130 

.128 

13.3 


» Volucs of vfacoslty of air obtained from linear variation with temperature given In ref. 16. 
Calculated as indicated In appendix D. 


The rate of total water catch on the airfoil per unit span is 
calculated from 

TF„=0.3296A^„T7iW 

where T is in feet, w in grams per cubic meter, and Vi in 
miles per hour. For this particular example, the table shows 
that for the 4-foot chord, the double-wedge airfoil has a 
somewhat larger value of E„ than the NACA 0006-64 air- 
foil; and for the 20-foot chord, the opposite is the case. 
Comparison of the two airfoils for a given chord length 
shows that, as expected, the effect of the total collection 
efficiency on the rate of water catch is small when compared 
with the free-stream velocity ratio chosen (2:1). The most 
significant comparison to be obtained from the table is that, 
for the double-wedge airfoil, increasing the chord length by a 
factor of 5 (from 4 to 20 ft) results in an increase of only 28 
percent in the rate of water catch; whereas, for the NACA 
airfoil, a like change in the chord length results in an increase 
of 94 percent. This difference is the result of the fact that, 
unlike the local impiogement on the rounded-leading-edge 
airfoils at free-stream Mach numbers less than critical, the 
local impiugement at a given point on the surface of a 
double-wedge airfoil does not vary with the chord length of 
the airfoE, which is at a supersonic Mach number above the 
shock-wave-attachment Mach number. 


APPENDIX D 


TOTAL 


COLLECTION EFFICIENCY OF 


SYMMETRICAL DOUBLE-WEDGE AIRFOILS AS FUNCTION OF MODIFIED INERTIA 
PARAMETER 


The total collection efficiency E„ is presented in figure 11 
as a function of the modified inertia parameter Kq, as sug- 
gested in reference 14. The relation for Ko is 

• (25) 


where K is the inertia parameter, defined as 


^_^2 

9 nic 


(24) 


Since K is equal to X,/c (ref. 4), may also be written as 



The free-stream droplet Reynolds numbers used herein in 
calculating range up to 2104. The empirical drag law 
(eq. (5)) used in this report for the droplet motion in the 
moving reference frame is vahd up to the Reynolds number of 
approximately 600. Therefore, for the values of Kfj pre- 
sented in figure 11 (and also fig. 10), the values of X/X, were 
obtained from table I of reference 4. 

The curves of E„ against Si, in figure 11 have exactly the 
same form as the curves of E„ against Fo, except that the use 
of Si) results in a displacement of the entire curve toward 
larger abscissa values. The displacement of the S„ against 
Si) curves obtained by varying the thickness ratio (fig. 1 1 (b)) 
or the free-stream Mach number (fig. 11 (c)) is less than 
the displacement of the E„ against Fq curves obtained by the 
same procedures (figs. 9 (b) and (c)). 
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